The Spin-SAF transition in NaV^Os induced by spin-pseudospin coupling 



o 
o 



Claudius Gros and Gennady Y. Chitov 
Fakultat 7, Theoretische Physik, University of the Saarland, 66041 Saarbriicken, Germany. 

(Dated: February 2, 2008) 

We present microscopic estimates for the spin-spin and spin-speudospin interactions of the quarter- 
filled ladder compound NaV2 0s, obtained by exactly diagonalizing appropriate clusters of the un- 
derlying generalized Hubbard Hamiltonian. We present evidence for a substantial interladder spin- 
pseudospin interaction term which would allow simultaneously for the superantiferroelectric (SAF) 
charge (pseudospin) ordering and spin dimerization. We discuss the values of the coupling constants 
appropriate for NaV20s and deduce the absence of a soft antiferroelectric mode. 

PACS numbers: 75.30.Gw,75.10. Jm,78.30.-j 
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Introduction- The insulating transition-metal com- 
pound NaV205 is, until now, the only known quarter- 
filled symmetric ladder compound the closely re- 
lated LiV205 has asymmetric ladders |2|. It has a sin- 
gle electron per rung which is located in the intra-rung 
V 45 — V 4 ' 5 bonding orbital, forming spin- 1/2 Heisenberg 
chains in the high-temperature state |3(. NaV^Os has a 
dual phase transition at T c = 34 K. Below this temper- 
ature two things happen: A spin gap opens |4|, like in 
a spin-Peierls system and a crystallographic lattice dis- 
tortion occurs which leads to a charge disproportionation 
y 4.5-8 _t/4.5+<5 on ^- ne run g Sj alternating along the ladder 
direction 0, 0> S 13 • Note also that in the literature on 
NaV205 its charge order is called the "zig-zag phase", 
what characterizes the antiferroelectric order in a sin- 
gle ladder only. In fact, the two-dimensional long-range 
charge order in NaV^Os is the super-antiferroelectric 
(SAF) 11], and we proposed to call the transition in 
NaV205 the spin-SAF transition. To develop an under- 
standing of this fascinating phase transition is an impor- 
tant matter, as it may serve as a model for other systems 
with coupled spin and orbital degrees of freedom|9j. 

Basing ourselves on the earlier suggestion of Mostovoy 
and Khomskiij^, l| that a bilinear coupling between the 
charge and spin degrees of freedom, similar to the spin- 
phonon coupling in spin-Peierls systems, may be respon- 
sible for the transition in NaV^Os, we have proposed the 
theory of the spin-SAF transition in the coupled spin- 
pseudospin which we believe can ex- 
plain the situation in that compound. In our theory the 
simultaneous appearance of the SAF charge charge order 
and of spin dimerization is driven by a single interlad- 
der spin-pseudospin coupling. Here we present results 
from a microscopic study for the strength of the coupling 
parameters necessary to understand the spin-SAF phase 
transition in NaV205. 

Microscopic Hamiltonian- The quarter-filled ladder 
compound NaV2 05 is described microscopically by a 
generalized Hubbard model |l2| with hopping matrix- 
elements t a and tb across the rung of the ladder (the crys- 
tallographic a-direction) and along the leg of the ladder 
respectively (the b-direction) , compare Fig. ^ in addi- 
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FIG. 1: The clusters used for the exact-diagonalization study, 
(a) A 2x4 ladder with four electrons with are dominantly lo- 
cated on the V-V bonding orbital, indicated symbolically by 
the dashed ellipses, (b) A six-site cluster with two rungs and 
two sites of the respective adjacent ladders, with four elec- 
trons, (c) Same as (b) but with two electrons. 



tion to an interladder hopping matrix element t±. The 
Coulomb interaction between the electrons gives rise to 
the onsite repulsion U, and the intraladder n.n. matrix 
elements V a (rung) and V& (leg). Here we neglect the 
possible intraladder diagonal repulsion Vd and further in- 
terladder matrix elements. 

Pseudospin Hamiltonian- NaV^Os is an insulator and 
the low-energy excitations of the microscopic generalized 
Hubbard model can therefore be mapped in perturbation 
theory 0,0], for small tb/U and small tb/V a , to a spin- 
pseudospin Hamiltonian. Here we consider it in the form 
H = H T + H s + H S t, with 



Ht 
Hst 



n,m m,n 

Jl 5 ^n,m ' ^n+l,m (2) 

n,m 

£ S n>m ■ S n+ i :?n (T n rn+1 — T n jfn _i) (3) 



where the spin/pseudospin operators obey the usual spin- 
algebra, e.g. [T a ,T 13 ] — iea^^T 1 . The sites indices 
(n, m) count rungs/ladders. In the high-temperature 
phase of NaV^Os the electrons are delocalised on the 
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FIG. 2: Estimates of the interladder spin-pseudospin coupling 
e, in units of the estimated intraladder spin-exchange Ji, for 
t/ = 3, to, = 0.35, it = 0.15 by comparing the respective 
triplet-excitation energies of the six-site clusters (b) and (c) 
in Fig. □ 



rungs with (T£ m ) ~ —1/2. in the low-temperature phase 
a finite-degree of charge ordering occurs with (T* rn ) cx 
(-!)"• 

In @ we have neglected further intraladder spin- 
orbital coupling terms 13] which are not critically rel- 
evant for the phy sics of the spin-SAF phase transition in 
NaV 2 5 [HEI. The 

magnitude of the spin-pseudospin 
coupling term e in 10 has not yet been estimated in per- 
turbation theory, it is ~ 

Recently we have shown [l(l , that the spin- pseudospin 
coupling term Hst ~ £ could at the same time drive 
the observed phase transition in NaV205 and lead to the 
observed opening of a gap in the spin-excitation spectrum 
via an alternation of the effective spin-spin coupling 
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along the ladder. 

Spin-pseudospin coupling- The exchange coupling J\ 
in can be estimated by diagonalizing small clusters 
of H and comparing the energy splitting between the 
ground-state singlet and the first excited triplet with the 
corresponding energy gap of the respective finite-size spin 
model. 

The inter-ladder spin-pseudospin coupling e can be de- 
termined using First we calculate the singlet-triplet 
energy gaps AE t (N) of the six-site clusters illustrated in 
Fig. (b) and (c) respectively i.e. with N — 4, 2 elec- 
trons. The position of the spins inside the clusters shown 
in Fig. Q illustrate typical spin-configurations realized 
for the parameters relevant for NaV^Os. AE t (N) there- 
fore correspond to the triplet-energy gaps of a system 
J cff Si • S2 with two spins Si and S2 located on two n.n. 
rungs of the ladder, which is just J cS . For the six-site 
cluster with N = 4,2 electrons we have approximatively 
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FIG. 3: Ratio A s /A t of the singlet (A s ) and the triplet 
(At) gap of spin-1/2 clusters with periodic boundary con- 
ditions and n.n./n.n.n. couplings Ji and J2 respectively. For 
J2/J1 > 0.24 this ratio extrapolates to zero in the thermody- 
namic limit. 



Using J|} we then find 

e = (A£ t (2)-A£ t (4))/2 (5) 

for the spin-pseudospin coupling e. The results are pre- 
sented in Fig. 10). 

Parameters for NaV^On- The Coulomb repulsion U 
has been estimated [12, Ej to be U w (3 - 4.1) eV for 
NaV2C>5. Its exact value is not critical and we use here 
U = 3.0 eV, all units throughout this papers are in eV. 
We use t a — 0.35 eV for the intra-rung hopping t a , close 
to the LDA-estimate H El of 0.38 eV. The inter-rung 
hopping if, along the leg of the ladder is 0, Q approx- 
imately 0.15 eV - 0.175 eV. 

The absolute magnitude of the interladder hopping 
matrix element t\ has been difficult to determine from 
first principles. It is relatively small @,H>E1 with a 
magnitude up to the estimated total interladder coupling 
of the order of 0.06 eV. Additionally, one has attributed 
the smallness of the frustrating inter-ladder spin-spin 
coupling in the high-temperature state 0, EH to a 
mutual cancellation of antifcrromagnetic and ferromag- 
netic contributions due to intermediate singlet and triplet 
states EI El- The substantial size of e found in our 
study and presented in Figgis then in qualitative agree- 
ment with the finding of Yaresko et al. of a substantial 
ferromagnetic interladder spin-spin coupling, within a 
spin-resolved density-functional study of the fully charge- 
ordered state We note that the true amount of 
charge disproportionation (T* m+1 — T* m _ 1 ) in NaV20s 
will be substantially reduced from unity and with a cor- 
responding reduction of the size of the intraladder dimcr- 
ization via Q. 

The size of the inter-size Coulomb repulsion matrix el- 
ements V a and Vb have been estimated to be up to half an 
eV EH i with a reduction from the respective bare values 
due to screening. The Ising-like pseudospin coupling- 
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FIG. 4: The frustration parameter a — J2/J1 for the 
n.n./n.n.n. couplings Ji and J2 along the ladder, estimated by 
comparing A s /A t for the 4-site Heisenberg ladder (see Fig. |HJ 
with the results obtained for the 2x4 ladder (see Fig. Q^ a )) 
for U = 3, t a — 0.35 and various V a and 14.. For NaV^Os the 
hopping tb along the leg is about tb ~ 0.175. 



constant g in is given by g = 4(V& — V4), where Vd 
the diagonal Coulomb repulsion matrix clement in be- 
tween two rungs 0. Due to geometry we have roughly 
V a ~ Vfc 1=3 V^Vd- In order to avoid a proliferation of 
parameters we have set here Vd — and examined two 
cases, see Fig. H = V a and V = 14/3. The latter case 
simulates the reduction of g » 4(14 — Vb/Vty ~ 414/3 
by V d . 

Intraladder frustration- Vojta et al. [17| have stud- 
ied recently the possibility of a intraladder spin gap for- 
mation for extended quarter filled Hubbard models on 
a single two-leg ladder, similar to the one studied here. 
They find a spontaneous spin gap formation by DMRG 
for tb > t a - This value for tb is outside of the param- 
eter range relevant for NaV20s but of general interest. 
Here we present a novel method to analyze finite-size data 
which allows to obtain directly the microscopic frustra- 
tion parameter a = J2/J1, where Ji is the n.n.n. spin- 
coupling. 

The frustrated spin-1/2 Heisenberg chain sponta- 
neously dimerizes for a > a c and the critical a c rs 0.24 
[Til can be accurately determined by the crossing of the 
singlet excitation gap AE S (L) with the triplet excitation 
gap AE t (L) in chains with length L. For a < a c we have 
AEt(L) < AE S (L), for a > a c the other way around. 

In Fig. H we plot the ratio AE S (L) / AE t (L) for some 
small Heisenberg chains of length L — 4, 8, 12, as a func- 
tion of J2/J1. The Majumdar-Gosh state is exactly real- 
ized for a = 0.5. For L — 4 there are only half as many 
J2-bonds than Ji-bonds and the Majumdar-Gosh state 
is realized for J2 = Ji , we have therefore plotted the data 
in Fig. |3]for L = 4 as a function of (J 2 /2)/Ji. 

The finite-size scaling properties of AE S (L) / AE t (L) 
are here, however, unimportant. Here we point out that 
the ratio AE S (L) / AE t (L) is characteristic for a given L. 



By calculating this ratio for the four-rung quarter-filled 
ladder illustrated in Fig. ^ we can therefor accurately 
determine the degree of intraladder frustration. The re- 
sults are presented in Fig. 01 We confirm that the system 
spontaneously dimerizes for % — t a — 0.35, for the pa- 
rameters V a — Vb considered by Vojta et al. |l7|]. We 
also find a tendency towards an anti-frustration, i.e. a 
ferromagnetic Ji for intermediate values of J2, indicat- 
ing the absence of frustration for NaA^Os. In the limit 
V a = = Vb we find a very large and negative J2 indi- 
cating long-range interactions. Here the analysis breaks 
down, as it is valid only for small to moderate values of 
Ji. 

Finally we present in Fig. El our estimates for the n.n. 
intraladder spin coupling Ji, obtained by comparing the 
triplet-gap of the four-rung ladder (see Fig. with that 
of the four-site Heisenberg chain. For % w 0.15 — 0.17, 
the range relevant for NaV20s we find very reasonable 
value J 1 = 0.04 — 0.06, close to the experimental value Q 
of 560 K = 0.048 eV, which are far more accurate than 
those by simple perturbation theory 

Critical temperature- From the theory of the spin-SAF 
transition in the coupled spin-pseudospin model [ToL lisj 
we obtain the following equation for the critical temper- 
ature at the Ising couplings g < g a (i.e., at the couplings 
where the pure transverse Ising model Q is always dis- 
ordered) 



viTc/Jx) = Ji (5- ffo)/4£ 2 



(6) 



where T](x) is the dimerization susceptibility of the 
dimerized Heisenberg XXX chain with the Hamiltonian 
ffxxx = E„ -Ml + (-l) n A)S n S n+1 in the limit A -» 0, 
and <?o = 4i a is the mean-field approximation for the 
quantum critical point of the transverse Ising model Q). 
The leading term of dimerization susceptibility, as fol- 
lows from the direct numerical (DMRG) calculations of 
Kliimper and co-workers 20J, can be reasonably approx- 
imated by the bosonization result of Cross and Fisher 
|U, i.e., rj(x) « a /x, where we take a Q ~ 0.26. 
Using the above approximation for r)(x), we obtain the 
critical temperature 



4a e 2 /(5o - .9) 



(7) 



as a function of the couplings. 

Parameters for NaV2 0s~ Taking the experimental 
T c = 34 K = 2.93 meV and a spin-pseudospin coupling 
of about £ w 0.4 Ji = 19.3 meV we find from Eq. 
that g - g w 0.132 eV, or g/g w 0.91. This es- 
timates for the intrachain Ising coupling implies that 
NaV2C>5 is on the disordered side, relatively close to the 
mean-field quantum critical point of the transverse Ising 
model. The standard RPA 0, |2^| for the transverse 
Ising model gives the following spectrum of the pseu- 
dospin (charge) excitations Eq — y / 2t a [2t a + m x g(q)], 
where the pseudospin average m x = tanh(t a /T)/2 w 
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FIG. 5: Estimates for the n.n. spin-coupling Ji along the 
ladder, see Eq. © from exact diagonalization of the 2x4 
ladder, see Fig.^a). The parameters are for U = 3, t a — 0.35 
and various values for V a and V&. For NaV2 0s the hopping 
tb along the leg is about tb ~ 0.175. 



1/2 at T <C i a and g(q) is the Fourier transform 
of the Ising coupling. _E q corresponds to the pole of 
the retarded pseudospin-pseudospin correlation function 
(T z T z )(co,q) = D zz (u,q) = 2t a m x /(cJ 2 - E\) @. In 
the function g(q) found earlier on the effective 2D lattice 
we retain only the largest intraladder coupling g, so 

(s < go) 



^o/qsAF = ztaV 1 ±g/go ■ 



(8) 



With the above estimates for the couplings we get Eq w 
0.97 eV » 7800 cm" 1 which agrees well with the ob- 
served peak of the optical conductivity (q rs 0) |2^| . At 
q = qsAF corresponding to the SAF ordering wave- vector 
(qsaf = (0, tt) on the effective 2D square lattice [l^) we 
get Eq SAF w 0.21 eV « 1700 cm" 1 . 

Discussion- From the Cross-Fisher RPA for the cou- 
pled spin-phonon system [2l| , Gros et al have shown 
that a soft mode (i.e., a pole in the phonon-phonon cor- 
relation function) occurs at a spin-Peicrls transition only 
if the energy of the involved phonon mode f2 is small, 
namely, < 2.2 T$p- The RPA for the coupled spin- 
pseudospin model renormalizes the pseudospin correla- 
tion function as D~ z R = D" 1 — II with II = —is 2 ?], 
resulting in Eq. 10 for the critical temperature of the 
spin-SAF transition and the same condition for the oc- 
currence of a pole of D ZZt R. A soft antiferroelectric mode 
is therefore absent in NaV 2 5 , since -E qsAF w 210 meV > 
2.2 T c = 6.4 meV. 

The energy 2i? qsAF » 3400 cm" 1 is quite close to the 
position of the broad peak near 4000 cm -1 observed in 
the optical absorption experiments at the room temper- 
ature (see also for more detailed data in the fre- 
quency range concerned). This suggests that this peak 
can be attributed to the two-particle (pseudospin) exci- 
tations with the ordering wave-vectors. 



Conclusions- We have proposed that the spin-SAF 
transition is of novel kind, in the sense that it is driven 
directly by the coupling in between spin and orbital de- 
grees of freedoms. We have presented support for this 
scenario from numerical estimates for the coupling pa- 
rameters and proposed an mechanism for a hitherto un- 
explained feature in the infrared spectrum. These stud- 
ies might be of relevance for other materials with strong 
spin-orbital couplings. 
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